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Finding control fields (pulse sequences) that can compensate for the dispersion in the parameters 
governing the evolution of a quantum system is an important problem in coherent spectroscopy and 
quantum information processing. The use of composite pulses for compensating for dispersion in 
system dynamics is widely known and applied. In this paper, we introduce new pulse elements for 
correcting pulse errors. These design methods are analytical and can be used to prove arbitrarily 
good robust performance. Furthermore, the time-to-compensation is superior to existing Fourier 
Synthesis Methods which is critical for minimizing errors due to relaxation effects. 



I. INTRODUCTION 

Many applications in control of quantum systems in- 
volve controlling a large ensemble by using the same con- 
trol field. In practice, the elements of the ensemble could 
show variation in the parameters that govern the dynam- 
ics of the system. For example, in magnetic resonance 
experiments, the spins of an ensemble may have large dis- 
persion in their natural frequencies (Larmor dispersion), 
strength of applied rf-field (rf-inhomogeneity) and the re- 
laxation rates of the spins [2] to name a few. In solid state 
NMR spectroscopy of powders, the random distribution 
of orientations of inter-nuclear vectors of coupled spins 
within an ensemble leads to a distribution of coupling 
strengths |13j . A canonical problem in control of quan- 
tum ensembles is to develop external excitations that can 
simultaneously steer the ensemble of systems with vari- 
ation in their internal parameters from an initial state 
to a desired final state. These are called compensating 
pulse sequences as they can compensate for the disper- 
sion in the system dynamics. From the standpoint of 
mathematical control theory, the challenge is to simul- 
taneously steer a continuum of systems between points 
of interest with the same control signal. Typical appli- 
cations are the design of excitation and inversion pulses 
in NMR spectroscopy in the presence of larmor disper- 
sion and rf-inhomogeneity [2HH1 ITUl - fT2"] or the transfer of 
coherence or polarization in coupled spin ensemble with 
variations in the coupling strengths [14] . In many cases 
of practical interest, one wants to find a control field that 
prepares the final state as some desired function of the 
parameter. For example, slice selective excitation and 
inversion pulses in magnetic resonance imaging [15H18) . 
The problem of designing excitations that can compen- 
sate for dispersion in the dynamics is a well studied sub- 
ject in NMR spectroscopy and extensive literature exists 
on the subject of composite pulses that correct for dis- 
persion in system dynamics 2-8 . Composite pulses have 
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recently been used in quantum information processing to 
correct for systematic errors in single and two qubit op- 
erations [TM2"3] . 

The focus of this paper is to present novel pulse el- 
ements for compensating errors arising from uncertain- 
ties/imperfections in the pulse amplitude. The construc- 
tions presented here have the advantage that they are 
analytical and exhibit favorable performance compared 
to the existing analytical Fourier Synthesis Method [T|. 
Namely, a higher level of rf robustness is obtained for 
the same pulse length and power, which is critical for 
minimizing errors due to relaxation. 

To fix ideas, consider an ensemble of noninteracting 
spin i particles in a static field Bo along the z axis and a 
transverse rf-field, (A(t) cos(0(i)), A(t) sin(0(i))), in the 
x — y plane. Let x, y, z represent the coordinates of the 
unit vector in the direction of the net magnetization vec- 
tor of the ensemble. The dispersion in the amplitude of 
the rf-field is given by a dispersion parameter e such that 
A(t) = eA (t) where e G [1 - S, 1 + S], for S > 0. The 
Bloch equations for the ensemble take the form 
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where 

(u(t),v(t)) = j(A (t) cos(0(t)), A (t) sin(0(<))). 

Let X(t) denote the unit vector (x(t),y(t), z(t)). Con- 
sider now the problem of designing controls u(t) and v(t) 
that simultaneously steer an ensemble of such systems 
with dispersion in the strength of their rf-field from an 
initial state X(0) = (0, 0, 1) to a final state X F = (1, 0, 0) 
|10j . This problem raises interesting questions about con- 
trollability, i.e., showing that in spite of bounds on the 
strength of the rf-field, yju 2 {t) + v 2 (t) < A max , there 
exist excitations (u(t),v(t)), which simultaneously steer 
all the systems with dispersion in e, to a ball of desired 
radius r around the final state (1,0,0) in a finite time 
(which may depend on A maxi B, S, and r). Besides steer- 
ing the ensemble between two points, we can ask for a 
control field that steers an initial distribution of the en- 
semble to a final distribution, i.e., different elements of 



2 



the ensemble now have different initial and final states 
depending on the value of the their dispersion parameter 
e. The initial and final state of the ensemble are described 
by functions -Xo(e) and Xp(e) respectively. Consider the 
problem of steering an initial distribution Xo(e) to within 
a desired distance r of a target function Xp{e) by appro- 
priate choice of controls in equation ([I]). We use the L2 
norm as our error metric 

rl+S 

E 2 = / \\X F (e)-X target (e)\\ 2 de. (2) 
Ji-s 

If a system with dispersion in its parameters can be 
steered between states that have dependency on the dis- 
persion parameter arbitrarily well, then we say that the 
system is ensemble controllable with respect to those pa- 
rameters. 

This paper is organized as follows. In the following sec- 
tion, we introduce the key ideas and through examples, 
highlight the role of non-commutativity in the design of 
a compensating control. We present novel pulse elements 
that compensate for inhomogeneity or uncertainty in the 
amplitude of the rf-field. The presented method extends 
known techniques for pulse sequences that are robust to 
rf inhomogeneity. The methods presented may also find 
applications in design of anisotropy compensating pulse 
design or in solid state NMR experiments [9]. 



To effect a larger amplitude rotation, we consider the 
sequence of transformations 

Hx = l[(V k ) nk ~ exp j e £ a k cos{ lk e)n y j , (7) 

k \ k I 

where a k = n[3 k . In practice, (3 k < jg is suitably small 
and results in an error that is less than 1% in the L2 
sense in Now, the coefficients a k can be chosen so 
that 

^2 a k cos(7 fc e) w - (8) 

k 

for 1 - 8 < e < 1 + S, with < S < 1. 
Therefore, 

H 1 ~ exp(6n y ) (9) 

approximately independent of e. The dependence on e 
can be made arbitrarily small by increasing the pulse 
length and extending the number of terms in the sum- 
mation leading to pulses that are immune to dispersions 
in the rf amplitude as claimed. 



B. Modified Fourier Synthesis Using S Modulation 



II. NOVEL PULSE ELEMENTS FOR 
COMPENSATING RF-INHOMOGENEITY 

A. Fourier Synthesis Technique 

Example 1: Main Concept 

To fix ideas, we begin by considering the Bloch equations 
in a rotating frame with only rf-inhomogeneity and no 
Larmor dispersion. 



X = e(u(t)Vt y +v(t)Q. x )X, 
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are the generators of rotation around x, y and z axis, 
respectively. 

We define the pulse elements pQ 

U\{Pk) = exp(-7fcfi x e) ex P( ^ ^V) exp^f^e), (4) 
U 2 {Pk) = exp(j k Q x e) exp(^-flye) exp(-^ k fl x e). (5) 

which correspond to directly accessible evolutions. 
For suitably small j3 k , we have 



V k = U 2 Ui ~ exp(e/3 fc fij,cos(7 fc e)). 



(G) 



In this section we develop a modified Fourier Synthe- 
sis Method that will be shown to have favorable time- 
robustness properties to the original Fourier Synthesis 
technique. To this end we consider the following system 



Y = (eu(t)n x +v(t)n z )Y 



(10) 



which corresponds to a system with one pure control 
by way of v(t) and one control with dispersion, eu(t). 
We will apply a similar Fourier Synthesis Method (FSM) 
analysis on the system and show that this results in a 
modified Hamiltonian to that of the previous section, the 
advantages of which will be discussed in section III. Wc 
then show how the previous system with both controls 
exhibiting dispersion ([3]) can be transformed into ( 10 1 by 
an appropriate change of coordinates. 
Consider the modified transformations 

Ui = exp(~j k efl x ) exp( y^ 2 ) exp(j k en x ) (11) 

U 2 =exp(7 fe en s )exp(--^fi z )exp(-7 fe efL c ) (12) 

By again choosing fi k sufficiently small, we have 

V k = U 2 Ui ~ exp(/3 fc sin(7 fc e)ft a ). (13) 
Applying a sequence of such transformations 



H 2 = Y[(v k ) nk 

k 



exp a k sin(7 fc e)f2 



v > 



(14) 
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FSM Expansion With Two Terms 



8 Mod Expansion Two Terms 




■ Target 7t/2 

■2 Term Expansion 



Epsilon (e) 
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FIG. 1. Left: 2 Term approximation for Hi = Right: 2 Term approximation for Hi = f using two terms with gradient 
descent for frequency selection 



where again ctk = nftk is used to control the error from 



the approximation in (13). Now, the coefficients and 



7/c can be chosen so that 



u k sin(7 fc e) 



(15) 



over the range of e of interest 1 — <5<e<l + <5 resulting 
in a robust rotation. We point out that ( 15 ) resembles 



([8j), but no longer contains an e factor external to the 
trigonometric argument and that cos has been replaced 
with sin. 

To see how (10) can be generated from (|3j), we return 
to (|| 

X = eA{cos (0i(t) + 4> 2 {t)) tl x + sin^ifX) + fc(t)) U y )X 



which corresponds to ( 10 ) once the appropriate identifi- 
cations are made. 

This means that Vk can be directly produced by im- 
plementing (f>i(t) as 0, 7r, 7r and over At time intervals 
such that A At = jk, and with —02 a delta pulse with 
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and — ^ at time At and 3 At respectively in 



^. As <j) 2 {^At) = 0, X{AAt) = Y{4At) and the lab 
frame and Y frame coincide after each pulse sequence, 
completing the 5 Modulated Pulse Design Method. 

The only remaining task is selection of the amplitude 
and frequency coefficients which we defer until the follow- 
ing section where we discuss the associated performance 
of the algorithms. 

In this paper we have focused on producing e-robust 
uniform rotations, however the reader should note that 
by appropriate selection of in Hi and H 2 , any rotation 
as a function of e can be produced. 
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and move into the frame 

r = ex p(-<fe(t)n z )x 

Y = [eA(t)( cosr>i(£)O x 



- sin 



(16) 

0i(t)n„ )-0 2 ojy(i7) 



Heuristic 

Error FSM 
Error 5 Mod 
Flip Z FSM 



n=2 n=3 n=4 
0.06831 0.06523 0.06473 
0.02012 0.00290 0.00044 
127.120 187.200 199.600 



Greedy 

Error FSM 
Error S Mod 
Flip Z FSM 



III. SIMULATIONS AND ERROR 
PERFORMANCE OF F. SYNTHESIS AND S 
MODULATION 



n=2 n=3 n=4 
0.04031 0.01506 0.00941 
0.02029 0.00422 0.00247 
130.497 179.062 229.101 



Gradient 

Error FSM 
Error 5 Mod 
Flip Z FSM 



n=2 n=3 n=4 
0.07339 0.01874 0.00423 
0.01940 0.00280 0.00044 
120.519 225.780 347.413 



Flip Z 6 Mod 115.230 172.001 216.138 Flip Z 6 Mod 110.952 165.178 216.177 Flip Z r5 Mod 113.341 170.134 216.137 

Table 1: Performance of Fourier Synthesis and S modulation for heuristic, greedy and gradient descent based 
frequency selection. In all cases, 5 modulation outperforms the Fourier Synthesis method in terms of L2 Error for a 
given pulse duration, which is important for minimizing relaxation effects. L2 error is calculated with respect to the 
desired final magnetization [1,0,0]'. 
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FIG. 2. Left: L2 Error with respect to the desired final magnetization [1,0,0]' for two term Hi and H2 as a function of the 
dispersion parameter e where frequencies were found using gradient descent. Right: Resulting X-magnetization for the 2 term 
sequence for both Fourier Synthesis and S Modulation. 8 modulation has a favorable magnetization profile while requiring a 
shorter pulse duration. 



The previous section reduced the problem of RF dis- 
persion compensation to parameter selection, 7fc and a^. 
Given the inhomogeneity parameter, e G [1 — <5o, 1 + #o]j 
we compute the error performance of synthesizing the 
effective Hamiltonians 



Hi = exp {^Yl a k cos ( £ 7fc)^^ 
H 2 = exp a k MnD^y) 



(18) 
(19) 



to approximate 



H — exp(6tt y 



We note that the {ot\} can be directly calculated given 
the {7^} as 



1+5 



a 1 =$-V; (f,g)= / f(e)g(e)de (20) 
(cos( 7 ie),cos(7]e)); V % = (cos^e), -> (21) 



and similarly the {a 2 } can be calculated as 

a 2 = ^V; (/, g) = /(e) ff (e)de (22) 

*y = (sin( 7 ? e ), sin( 7 2 e)) ; V { = (sin( 7 ?e), 8) (23) 

so that the problem reduces to selecting the optimal fre- 
quencies. 

We report performance for three frequency selection 
methods, heuristically, greedy selection and gradient de- 
scent and show that 5 modulation outperforms Fourier 



Synthesis Methods for all frequency selection methods. 
Unless stated otherwise, the notion of optimal is with re- 
spect to L2 error for a given pulse duration. L2 error is 
calculated with respect to the desired final magnetization 
[1,0,0] and pulse duration is reported in total flip angle. 

As a starting point, we consider the problem of se- 
lecting the optimal amplitudes given known frequencies 
which we will choose heuristically. As sine obtains its 
maximum at 7r/2 and is relatively horizontal about this 
point, a natural selection for the frequencies in (151 is 

Ik 



(2fc+l) 



. Similarly, selecting 7^ in (8) to maximize 
flatness about e = 1 corresponds to 



ecos(7 fe e) 



cos(7 fc ) - j k sin(7 fe ) = 0. (24) 



e=l 



Numerically solving gives the first several 7^ = 
[.860,3.426,6.437]. 

The amplitude coefficients a were then calculated ac- 
cording to (20)- (23 1. The comparative performance of 



standard FSM, Hi, to S modulation, H 2 , is tabulated in 
table 1 and a complete description of the pulses is given 
in the appendix. 

An alternative algorithm is to sequentially select the 
frequencies employing a greedy algorithm, in which al- 
ready determined frequencies are held fixed, and only 
the newest frequency is optimized over. Explicitly we 

sequentially minimize the cost functions with respect to 

1/2 
Ik 



rl+S 

*i(7fe»-,7i) = / y^afeCOs(e7^) 
J 1-8 



de (25) 



^ 2 (7 fe 2 ,-,7i 2 ) 



rl+6 

/ ^a fc sin(e 7fc 2 )-^ de (26) 
J 1-8 11 
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again using (20 1- (23 1 to calculate the a. This was done 
using gradient descent and numerically calculating the 
necessary derivatives. Table 1 shows that the 5 modula- 
tion outperforms standard FSM. 

The most general method we applied (and best per- 
forming) was simultaneously optimizing F± and F% with 
respect to all frequencies using gradient descent, where 
derivatives were again calculated numerically. As with all 
descent schemes, there is concern that one merely obtains 
a local minima. Moreover, the problem of unspecified fre- 
quencies is how to project onto an overrepresented sub- 
space which is known to have local minima. To combat 
such issues we chose the optimal result after numerous 
starting points and note that the performance exceeds 
the other methods and the results are displayed in table 
1. 

As an example we consider the resulting parameters 
from optimizing a two term S modulated pulse using gra- 
dient descent 
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[88.6°, 265. l c 



[105.5°, 16.6°]. 



These are converted into a pulse sequence by hrst divid- 
ing large amplitudes of a 2 , into repeated sequences with 
smaller amplitudes according to (14) using a threshold 
value of 9° , which yields the modified parameters 

7 2 ' = [88.6°,. ..,88.6°, 265.1°,.. .,265.1°] 



12 times 



12 times 



2 times 



.4°,.. .,8.4° 



2 times 



Pulses are calculated as described in section II. B where 
pulse elements are 

[(7fc)o(27fc)iso°-a fc /2(7fc)o] 

with numbers inside the parentheses representing the flip 
angle and the subscripts, the phase. Applying to the 
modified parameters above yields the pulse sequence 

[(88.6)o(177.1) 175 .6(88.6) ] x12 [(265.1) (530.1) 175 . 9 (265.1) 

The performance is displayed in figure 2. The more terms 
kept in the series, the longer the sequences and overall 
pulse, but the more robust. 



IV. GENERAL MODULATION SCHEMES 

In many ways S modulation is the most natural choice 
as it has a nice correspondence with existing FSM's. 
However, other modulation schemes are possible and 
their analysis is warranted for the sake of completeness 
or in the event that abrupt phase adjustments in the RF- 
fields are not available. We begin by considering linear 
modulation. 
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FIG. 3. Trajectory in the interaction frame for linearly mod- 
ulated controls. The trajectory is the black path from 1-2 and 
back to 1, followed by the blue path from 3-4 and back to 3. 



Linear Modulation 



Returning to equation ( 10 \ with At — \ and 



< t < At 
vr At < t < 3At 
3Ai < t < 4Ai 



(27) 



(28) 



x2 



Y = (Aecosfa(t)Cl x - fafl z )Y 
Moving into the interaction frame 

Z = exp(- J a A cos fa (t)dt fl x )Y 
Z = — fa (cose J Acosfa(t)dt fl z 

+ sine / Acos fa(t)dt n y ) Z (30) 



(29) 



We note that under the assumptions 



cos fa(t)dt= I fa(t)dt = {) (31) 



the Z frame will agree with the Y frame which will agree 
with the lab frame at time T, so that it is sufficient to an- 
alyze the system in the interaction frame. Letting fa{t) 
be a linear modulation of the form 



—B < t < 2At 
B 2At < t < 4At 



(32) 
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FIG. 4. L2 Error (Left) and Total Flip Angle (Right) for Heuristic Frequency Selection for Fourier Synthesis and 8 Modulation. 
Error for a Fixed Pulse Duration is Smaller for S Modulation. 



Error for Greedy Frequency Selection 



Total Flip Angle for Greedy Frequency Selection 
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FIG. 5. L2 Error (Left) and Total Flip Angle (Right) for Greedy Frequency Selection for Fourier Synthesis and S Modulation. 
Both Duration and Error are Smaller for 5 Modulation. 



we can analyze the resulting rotation with the Peano- 
Baker Series 



$ = I 



= I 



H(t)dt+ / H(t) / H{a 1 )da 1 dt- 
ii Jo Jo 

2 



H(t)dt + o 



I + AB I sm(Aet)dm y + O 
Jo 

I+~(l-cos(7re)) + o(j 



= I+ 8 -§(l-S) + o(Sf+o(^ 



(33) 



which, to first order, has resulted in an evolution with the 
dispersion term reversed. Combining with a directly ac- 
cessible evolution of 8B ^ +l5 - > £l y will produce a pulse that 
is robust to first order in 5. Figure 3 displays the trajec- 
tory in the interaction frame for the linearly modulated 
pulse with e > 1 and provides the intuition for why the 
dispersion term is negated. 



Arbitrary Modulation Schemes 

Other modulation functions are also possible. Let 
\f(t)\ < B be such a candidate modulation, then choos- 
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FIG. 6. L2 Error (Left) and Total Flip Angle (Right) for Gradient Descent Frequency Selection for Fourier Synthesis and 8 
Modulation. Both Duration and Error are Smaller for 8 Modulation. 



ing which can be used to produce new dispersion dependen- 

cies and thereby robust pulses as was done in the linear 



t>2 = 



( f(t) < t < At 

I f(2At -t) At <t< 2At 



case. 



2At<t<3At (34) V. CONCLUSION 



-/(* - 2At) 

I -/(4Af - t) 3At<t< AAt 

We have presented a new method for pulse design in 
the presence of RF-inhomogeniety that extends existing 
will produce a rotation Fourier Synthesis Methods. The method displays supe- 

rior time-error properties to conventional Fourier Synthe- 
7L / B\ 2 S ' S ^ e th° c ' s ' This method is analytical and can be used 

I - 4 J A f(t) sm(Aet)dm y + O (—) (35) to produce arbitrarily robust performance. 
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VI. PULSE PARAMETER APPENDIX 



Heuristic FSM, n=2 



a = [187.3,33.8] 
7 = [49.3, 196.5] 

Pulse: [(49.3)o(4.5) 9 o(98.5)i 8 o(4.5) 9 o(49.3)o] x21 [(196.5)o(4.2) 9 o(393.0)i 80 (4.2) 9 o(196.5)o] x4 

Heuristic FSM, n=3 

a = [201.1,49.2,7.3] 
7 = [49.3,196.5,369.0] 
Pulse : [(49.3)o(4.4) 90 (98.5) 18 o(4.4) 90 (49.3)o] x23 

[(196.5)o(4.1) 90 (393.0) 18 o(4.1) 90 (196.5)o] x6 [(369.0)o(3.6) 90 (738.0) 18 o(3.6) 90 (369.0)o] xl 

Heuristic FSM, n=4 

a = [175.2903, 18.3977, -10.8059, -5.67454] 
7 = [49.3,196.5,369.0,546.0] 
Pulse : [(49.3)o(4.4) 9 o(98.5)i 80 (4.4) 9 o(49.3)o] x20 [(196.5)o(3.1) 9 o(393.0) 18 o(3.1) 9 o(196.5)o] x3 

[(369.0)o(-2.7) 9 o(738.0)i 8 o(-2.7) 9 o(369.0)o] x2 [(546.0)o(-2.8) 9 o(1092.1) 18 o(-2.8) 90 (546.0)o] xl 

Heuristic 6 Mod, n=2 

a = [105.5,16.7] 
7 = [90, 270] 

Pulse: [(90.0)o(180.0)i 7 5.6(90.0)o] xl2 [(270.0)o(540.0)i 75 . 8 (270.0)o] x2 

Heuristic 6 Mod, n=3 

a = [108.3,22.4,4.3] 
7 = [90,270,450] 

Pulse : [(90.0)o(180.0)i 7 5. 8 (90.0)o] xl3 [(270.0)o(540.0)i 76 . 3 (270.0)o] x3 [(450.0)o(900.0)i 77 . 9 (450.0)o] xl 
Heuristic 8 Mod, n=4 

a = [109.8,25.7,7.1,1.2] 
7 = [90, 270, 450, 630] 
Pulse: [(90.0)o(180.0) 17 5. 8 (90.0)o] xl3 [(270.0)o(540.0) 175 . 7 (270.0)o] x3 

[(450.0)o(900.0)i 76 . 4 (450.0)o] xl [(630.0)o(1260.0)i 79 .4(630.0)o] xl 

Greedy FSM, n=2 

a = [191.9,35.9] 
7 = [49.9, 192.7] 

Pulse: [(49.9) (4.4) 90 (99.9) 180 (4.4) 90 (49.9)o] x22 [(192.7) (4.5) 9 o(385.4) 180 (4.5) 90 (192.7)o] x4 

Greedy FSM, n=3 

a = [197.4,40.9,-3.8] 
7= [49.9,192.7,502.9] 

Pulse: [(49.9)o(4.5) 9 o(99.9) 18 o(4.5) 9 o(49.9)o] x22 [(192.7)o(4.1) 9 o(385.4) 180 (4.1) 9 o(192.7)o] x5 
[(502.9)o(-1.9) 9 o(1005.8)i 8 o(-1.9) 90 (502.9)o] xl 



Greedy FSM, n=4 

a = [200.7,43.7,-5.9,-1.9] 
7= [49.9,192.7,502.9,666.8] 
Pulse: [(49.9)o(4.4) 9 o(99.9) 18 o(4.4) 9 o(49.9)o] x23 [(192.7)o(4.4) 90 (385.4) 18 o(4.4) 90 (192.7)o] x5 

[(502.9)o(-3.0) 9 o(1005.8)i 8 o(-3.0) 9 o(502.9)o] xl [(666.8)o(-0.9) 9 o(1333.7)i 8 o(-0.9) 90 (666.8)o] xl 

Greedy 5 Mod, n=2 

a = [105.5,16.6] 
7 = [86.7,259.1] 

Pulse: [(86.7)o(173.4) 175 . 6 (86.7)o] xl2 [(259.1)o(5 1 8.1) 17 5.8(259.1)o] x2 

Greedy 5 Mod, n=3 

a = [108.2,22.2,4.1] 
7 = [86.7,259.1,427.8] 

Pulse : [(86.7)o(173.4) 17 5. 8 (86.7)o] xl3 [(259.1)o(518.1)i 76 . 3 (259.1)o] x3 [(427.8)o(855.7)i 77 . 9 (427.8)o] xl 
Greedy 6 Mod, n=4 

a = [108.5, 22.9, 4.6, -.3] 
7 = [86.7,259.1,427.8,730.2] 
Pulse: [(86.7)o(173.4) 175 . 8 (86.7)o] xl3 [(259.1)o(518.1) 176 . 2 (259.1) ] x3 

[(427.8) (855.7)i 77 . 7 (427.8)o] xl [(730.2)o(1460.5)i 8 o. 2 (730.2)o] xl 

Gradient Descent FSM, n=2 

a= [163.4,-15.7] 
7 = [51.5,373.7] 

Pulse : [(51.5)o(4.3) 9 o(103.0)i 8 o(4.3) 9 o(51.5)o] xl9 [(373.7)o(-3.9) 90 (747.4) 18 o(-3.9) 9 o(373.7)o] x2 

Gradient Descent FSM, n=3 

a = [169.6,-23.9,-10.3] 
7 = [52.4,379.1,550.3] 

Pulse : [(52.4)o(4.5) 9 o(104.9)i 8 o(4.5) 9 o(52.4)o] xl9 [(379.1)o(-4.0) 90 (758.2) 18 o(-4.0) 9 o(379.1)o] x3 
[(550.3)o(-2.6) 9 o(1100.6)i 8 o(-2.6) 90 (550.3)o] x2 

Gradient Descent FSM, n=4 

a= [174.4,-30.6,-19.0,-5.1] 
7 = [53.1,381.4,554.0,727.9] 
Pulse: [(53.1)o(4.4) 9 o(106.1) 18 o(4.4) 90 (53.1)o] x20 [(381.4) (-3.8) 9 o(762.7) 18 o(-3.8) 9 o(381.4)o] x4 

[(554.0)o(-3.2) 90 (1108.0)i 8 o(-3.2) 90 (554.0)o] x3 [(727.9)o(-2.6) 9 o(1455.8)i 8 o(-2.6) 9 o(727.9)o] xl 

Gradient Descent 5 Mod, n=2 

a = [105.5,16.6] 
7 = [88.6,265.1] 

Pulse: [(88.6) (177.1) 17 5. 6 (88.6)o] xl2 [(265.1)o(530.1)i 75 . 9 (265.1)o] x2 

Gradient Descent 5 Mod, n=3 

a = [108.3,22.4,4.3] 
7 = [89.1,267.0,444.5] 

Pulse : [(89.1)o(178.1)i 7 5. 8 (89.1)o] xl3 [(267.0)o(534.1) 176 . 3 (267.0)o] x3 [(444.5)o(889.0)i 77 . 9 (444.5)o] xl 



Gradient Descent S Mod, n=4 



a = [109.8,25.7,7.1,1.2] 
7 = [90.0,270.0,450.0,630] 
Pulse: [(90.0) (180.0) 17 5. 8 (90.0) ] xl3 [(270.0)o(540.0) 175 . 7 (270.0)o] x3 

[(450.0)o(900.0)i 76 . 4 (450.0)o] xl [(630.0)o(1260.0)i 79 .4(630.0)o] 



